Bose-Einstein condensation in quasi2D trapped gases 
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We discuss BEC in (quasi) 2D trapped gases and find that well below the transition temperature 
Tc the equilibrium state is a true condensate, whereas at intermediate temperatures T < Tc one 
has a quasicondensate (condensate with fiuctuating phase) . The mean-field interaction in a quasi2D 
gas is sensitive to the frequency loq of the (tight) confinement in the "frozen" direction, and one 
can switch the sign of the interaction by changing ujo- Variation of loq can also reduce the rates of 
inelastic processes. This offers promising prospects for tunable BEC in trapped quasi2D gases. 
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The influence of dimensionality of the system of bosons 
on the presence and character of Bose-Einstein conden- 
sation (BEC) and superfluid phase transition has been a 
subject of extensive studies in spatially homogeneous sys- 
tems. In 2D a true condensate can only exist at T = 0, 
and its absence at finite temperatures follows from the 
Bogolyubov fc~^ theorem and originates from long-wave 
fluctuations of the phase (see, e.g., However, as 

was first pointed out by Kane and Kadanoff Q and then 
proved by Berezinskii Q , there is a superfluid phase tran- 
sition at sufficiently low T. Kosterlitz and Thouless ||] 
found that this transition is associated with the forma- 
tion of bound pairs of vortices below the critical tem- 
perature Tkt = {Trh^ /2m)ns (m is the atom mass, and 
Us the superfluid density just below Tkt)- Earlier the- 
oretical studies of 2D systems have been reviewed in |^] 
and have led to the conclusion that below the Kosterlitz- 
Thouless Transition (KTT) temperature the Bose liquid 
(gas) is characterized by the presence of a quasiconden- 
sate, that is a condensate with fluctuating phase (see [^). 
In this case the system can be divided into blocks with 
a characteristic size greatly exceeding the healing length 
but smaller than the radius of phase fluctuations. Then, 
there is a true condensate in each block but the phases 
of different blocks are not correlated with each other. 

The KTT has been observed in monolayers of liquid 
helium |Q. The only dilute atomic system studied thus 
far was a 2D gas of spin-polarized atomic hydrogen on 
liquid- helium surface (see for review). Recently, the 
observation of KTT in this system has been reported . 

BEC in trapped 2D gases is expected to be qualita- 
tively different. The trapping potential introduces a fi- 
nite size of the sample, which sets a lower bound for the 
momentum of excitations and reduces the phase fiuctua- 
tions. Moreover, for an ideal 2D Bose gas in a harmonic 
potential Bagnato and Kleppner pO[ | found a macroscopic 
occupation of the ground state of the trap (ordinary 
BEC) at temperatures T < « N^^'^huj, where N is 
the number of particles, and u the trap frequency. Thus, 



there is a question of whether an interacting trapped 2D 
gas supports the ordinary BEC or the KTT type of a 
cross-over to the BEC regime However, the criti- 

cal temperature will be always comparable with Tc of an 
ideal gas: On approaching Tc from above, the gas density 
IS Uc N/R^^, where i?T, « ^jTc/muo'^ is the thermal 
size of the cloud, and hence the KTT temperature is 
~ h^nc/m ~ N^/^huj « Tc. 

The discovery of BEC in trapped alkali-atom clouds 
|l2| stimulated a progress in optical cooling and trapping 
of atoms. Present facilities allow one to tightly confine 
the motion of trapped particles in one direction and to 
create a (quasi)2D gas. In other words, kinematically the 
gas is 2D, and in the "frozen" direction the particles un- 
dergo zero point oscillations. This requires the frequency 
of the tight confinement lo^ to be much larger than the 
gas temperature T and the mean-field interparticle in- 
teraction n^g (no is the gas density, and g the coupling 
constant). Recent experiments jl^-p^ indicate a realistic 
possibility of creating quasi2D trapped gases and achiev- 
ing the regime of quantum degeneracy in these systems. 
The character of BEC will be similar to that in purely 
2D trapped gases, and the main difference is related to 
the sign and value of the coupling constant g. 

In this Letter we discuss BEC in quasi2D trapped gases 
and arrive at two key conclusions. First, well below Tc the 
phase fluctuations are small, and the equilibrium state is 
a true condensate. At intermediate temperatures T < Tc 
the phase fluctuates on a distance scale smaller than the 
Thomas-Fermi size of the gas, and one has a quasicon- 
densate (condensate with fluctuating phase). Second, in 
quasi2D the coupling constant g is sensitive to the fre- 
quency of the tight confinement ujq and, for a negative 
3D scattering length a, one can switch the mean-field in- 
teraction from attractive to repulsive by increasing loq. 
Variation of luq can also reduce the rates of inelastic pro- 
cesses. These findings are promising for tunable BEC. 

In a weakly interacting Bose-condensed gas the corre- 
lation (heahng) length h/^mnog {g > 0) should greatly 
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exceed the mean interparticle separation. In (quasi) 2D 
the latter is ~ l/v^^Tmo, and we obtain a small parame- 
ter of the theory, {mg/2TTh^) <^ 1 (see [||). 

We first analyze the character of BEC in a harmoni- 
cally trapped 2D gas with repulsive interparticle interac- 
tion, relying on the calculation of the one-particle density 
matrix. Similarly to the spatially homogeneous case , 
at sufficiently low temperatures only phase fluctuations 
are relevant. Then the field operator can be written as 
^(R)=n|^(R)exp{i(/)(R)}, where 0(R) is the operator of the 
phase fluctuations, and no(R) the condensate density at 
T=0. The one-particle density matrix takes the form 



(*t(R)v[,(0)) = V«o(R)no(0)exp{-((5^(R))2)/2}. (I) 

Here 5(/)(R) = (/)(R) — 0(0), and R = at the trap center. 
For a trapped gas the operator (/)(R) is given by 



(R)=^[4no(R)] 



-1/2 r + 
J V 



h.c. 



(2) 



were is the annihilation operator of an elementary 
excitation with energy ty, and — ± Vi^ are the 
Bogolyubov u, v functions of the excitations. 

In the Thomas- Fermi regime the density no(R) has the 
well-known parabolic shape, with the maximum value 
n-Om = no{0) « {Nm/ngY^'^iL!, and the radius Rtf ~ 
(2/^/tow^)^/^. The chemical potential is = n^mg, and 
the ratio Tc//i ~ (jxTi' ImgYl"^ ^ 1. For calculating the 
mean square fluctuations of the phase, we explicitly found 
the (discrete) spectrum and wavefunctions of excitations 
with energies <C /i by using the method developed for 
3D trapped condensates ||l^ . For excitations with higher 
energies we used the WKB approach. At distances R 
greatly exceeding the wavelength At of thermal excita- 
tions (e^«r) near the trap center, for T^fi we obtain 
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■log {R/Xt). 
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We also find that Eq.(||) holds at any T for a homo- 
geneous gas of density nom, where at T ^ /i it repro- 
duces the well-known result (see Q). In a trapped gas 
for r <C /i, due to the contribution of low-energy exci- 
tations, Eq.(^) acquires a numerical coefficient ranging 
from 1 at i? ^ i?^ f to approximately 3 at i? ~ Rt f ■ 

The character of the Bose-condensed state is deter- 
mined by the phase fiuctuations at i? ^ Rtf- With 
logarithmic accuracy, from Eq.(H) we find 



((Wtf))^)«(|) (^V^'logiV. 



(4) 



In quasi2D trapped alkali gases one can expect a value 
^ I0~^ or larger for the small parameter mg/2'Kh^ ^ and 
the number of trapped atoms N ranging from 10* to 10^. 

Then, from Eq.(P we identify two BEC regimes. At 
temperatures well below the phase fluctuations are 



small, and there is a true condensate. For intermediate 
temperatures T < the phase fluctuations are large 
and, as the density fluctuations are suppressed, one has 
a quasicondensate (condensate with fluctuating phase). 

The characteristic radius of the phase fluctuations R^~ 
At exp (7r?i^/mT), following from Eq.(||) under the con- 
dition ((5(/)(R))^)^l, greatly exceeds the healing length. 
Therefore, the quasicondensate has the same Thomas- 
Fermi density profile as the true condensate. Correlation 
properties at distances smaller than R^ and, in particu- 
lar, local density correlators are also the same. Hence, 
one expects the same reduction of inelastic decay rates 
as in 3D condensates [^. However, the phase coherence 
properties of a quasicondensate are drastically different. 
For example, in the MIT type of experiment on inter- 
ference of two independently prepared quasicondensates 
the interference fringes will be essentially smeared out. 

We now calculate the mean-field interparticle interac- 
tion in a quasi2D Bose-condensed gas, relying on the 
binary approximation. The coupling constant g is in- 
fiuenced by the trapping field in the direction z of the 
tight confinement. For a harmonic tight confinement, 
the motion of two atoms interacting with each other via 
the potential V{r) can be still separated into their rela- 
tive and center of mass motion. The former is governed 
by V{r) together with the potential Vh{z) = mw^z^ 
originating from the tight harmonic confinement. Then, 
similarly to the 3D case (see, e.g. (l^), to zero order in 
perturbation theory the coupling constant is equal to the 
vertex of interparticle interaction in vacuum at zero mo- 
menta and frequency E = 2/x. For low E > Q this vertex 
coincides with the amplitude of scattering at energy E 
and, hence, is given by 



g = f{E) = j dvi;(v)V{r)rf{v). 



(5) 



The wavefunction of the relative motion of a pair of 
atoms, "0(1"), satisfies the Schrodinger equation 



-^ + V{v) + Vh{z) 



hujo 



il)[v) = Ei;(v). (6) 



The wavefunction of the free x,y motion ipfir) = 
ipo{z) exp(iqp), with Lpo{z) being the ground state wave- 
function for the potential Vff(z), p= {x,y}, and q = 
{2mE/h^y/^ . As the vertex of interaction is an analyt- 
ical function of E, the coupling constant for /i < is 
obtained by analytical continuation of f{E) to E < 0. 

The possibility to omit higher orders in perturbation 
theory requires the above criterion {mg /2'Kh?) <^ 1. 

The solution of the quasi2D scattering problem from 
Eq.(^ contains two distance scales: the extension of ip{r) 
in the z direction, I = (?i/mwo)^^^, and the characteristic 
radius i?e of the potential V{r). In alkalis it ranges from 
20 A for Li to 100 A for Cs. At low energies {qRe ^1) 
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the amplitude f{E) is determined by the scattering of 
the s-wave for the motion in the x, y plane. 

We first consider the limiting case Re- Then the 
relative motion of atoms in the region of interatomic in- 
teraction is not influenced by the tight confinement, and 
ip{r) in Eq.(||) differs only by a normalization cocfhcicnt 
from the 3D wavefunction: 



(7) 



In the interval, where i?e^f <C I, Eq.(^ takes the form 
ip = ''Pasir) = r/(^o(0)(l — a/'')- This expression serves as 
a boundary condition at r ^ for the solution of Eq.(H) 
with V{r) ~ {r ^ R^). The latter can be expressed 
through the Green function G(r,r') of this equation: 



ip{r) = (fioiz) exp(iqp) + AG{r, 0), 



(8) 



The coefficients A and rj are obtained by matching the 
solution (||) at r — > with ^a,s{i~)- 

Similarly to the case of a purely ID harmonic oscillator 
(see, e.g., [|o|), we have 

G(r 0)^- ^f exp{^(z^cott/4^2-g^^^t-V2+p74tP)} 
I Jo t^(47ri)3 sint 

Under the condition {^'t^hujo) at r^/ we obtain 



G 



1 



2(27r)3/2Z 



log 



(9) 



Omitting the imaginary part of G (^) and comparing 
Eq.(H) with -i/jas, we immediately find 



log 



(10) 



In Eq.(|) 



put V/ = <^o(0) = (l/27r/2)i/4_ xhen. 



using the well-known result / dr'03£)(r)y(r) = AttH a/ 



Eqs. 



and (|T^) lead to the coupling constant 
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m l/a+{l/V2n)log{l/Trq^P) 



(11) 



For fi < 0, analytical continuation of Eqs. (^ and ( |ll| ) 
to E = h^q-^ /m < leads to the replacement E ^ \E\ = 
2\fi\ in the definition of q. 

The coupling constant in quasi2D depends on q = 
(2TO|/i|/?i^)^/2 and, hence, on the condensate density. In 
the limit Z ^ |a| the logarithmic term in Eq.(pl]) is not 
important, and g becomes density independent. In this 
case the quasi2D gas can be treated as a 3D condensate 
with the density profile oc exp (— z^/Z^) in the z direction. 

As follows from Eq.(pT|), for repulsive mean-field inter- 
action in 3D (a > 0) the interaction in quasi2D is also 
repulsive. For a < the dependence g{l) has a resonance 
character (cf. Fig.l): The coupling constant changes sign 
from negative (attraction) at very large I to positive for 



I < = (|a|/\/27r) log (I/ttij^/^). This should describe 
the case of Cs, where a < —600 A [^,0 and the con- 
dition I ^ Re assumed in Eq. (pi]) is satisfied at I < h. 
Near the resonance point /, the quantity {m\g\/27rh?) be- 
comes large, which violates the perturbation theory for a 
Bose-condensed gas and makes Eqs. (||) and (^l|) invalid. 

For I < Re (except for very small I) we used directly 
Eqs. (^ and (^ and calculated numerically the coupling 
constant g for Li, Na, Rb, and Cs. The potential V{r) 
was modeled by the Van der Waals tail, with a hard 
core at a distance Rq <C Re selected to support many 
bound states and reproduce the scattering length a. The 
numerical results differ slightly from the predictions of 
Eq. (ITTl) . For Rb and Cs both are presented in Fig.l. 
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FIG. 1. The parameter mgj^-Kf? versus IjR^ at fixed no for 
Rb (a) and Cs (b). Solid curves correspond to the numerical 
results, dashed curves to Eq.(pj]). The dotted curve in (b) 
shows the result of Eq.(pj]) in the region where m\g\l2-Kfi^ ~ 1. 

The nature of the g(V) dependence in quasi2D can be 
understood just relying on the values of g in the purely 
2D and 3D cases. In 2D at low energies the mean-field 
interaction is always repulsive. This striking difference 
from the 3D case can be found from the solution of the 
2D scattering problem in |jl^ and originates from the 
2D kinematics: At distances, where Re ^ P ^ 
{q 0), the solution of the Schrodinger equation for 
the (free) relative motion in a pair of atoms reads ip oc 
\og{p/d)/ \og{l/qd) {d > 0). We always have increas- 
ing with p, unless we touch resonances corresponding to 
the presence of a bound state with zero energy (d — > oo). 
This means that it is favorable for particles to be at larger 
p, i.e. they repel each other. 

In quasi2D for very large I the sign of the interparticle 
interaction is the same as in 3D. With decreasing I, the 
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2D features in the relative motion of atoms become pro- 
nounced, which is described by the logarithmic term in 
Eg. (pi]) . Hence, for a > the interaction remains repul- 
sive, whereas for a < the attraction turns to repulsion. 

The obtained results are promising for tunable BEC in 
quasi2D gases, based on variations of the tight confine- 
ment and, hence, I. However, as in the MIT studies of 
tunable 3D BEC by using Feshbach resonances |^^, an 
"underwater stone" concerns inelastic losses: Variation 
of I can change the rates of inelastic processes. For opti- 
cally trapped atoms in the lowest Zeeman state the most 
important decay process is 3-body recombination. 

This process occurs at interparticle distances r ^ 
max{i?e,|a|} p3| , 3 . We will restrict ourselves to the 
case where Z ^ |a| and is also significantly larger than Re- 
Then the character of recombination collisions remains 3- 
dimensional, and one can treat them in a similar way as in 
a 3D gas with the density profile {no/y/nl) exp 

However, the normalization coefficient of the wavefunc- 
tion in the incoming channel will be influenced by the 
tight confinement. Relying on the Jastrow approxima- 
tion, we write this wavefunction as a product of the three 
wavefunctions •ip{rik), each of them being a solution of 
the binary collision problem Eq.(^). In our limiting case 
the solution is given by Eq.(^ divided by i^o(O) to recon- 
struct the density profile in the z direction. The outgoing 
wavefunction remains the same as in 3D, since one has a 
molecule and an atom with very large kinetic energies. 

Thus, in the Jastrow approach we have an additional 
factor r]'^ for the amplitude and r]^ for the probability 
of recombination in a quasi2D gas compared to the 3D 
case. Averaging over the density profile in the z direction, 
we can relate the quasi2D rate constant a to the rate 
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constant in 3D (see |24 for a table of a^u in alkalis): 



a = {if /TTf)a'iD- 



(12) 



As 77 is given by Eq. (|lo|) , for a > the dependence a(/) 
is smooth. For a < the rate constant a peaks at Z « 
and decreases as /"*/(/* — Z)^ at smaller I. This indicates 
a possibility to reduce recombination losses while main- 
taining a repulsive mean- field interaction (5 > 0). For 
Cs already at / « 200 A {U ~ 500A) we have a ~ 10"" 
cm'^/s, and at densities 10^ cm~^ the life-time r>l s. 

The predicted possibility to modify the mean-field in- 
teraction and reduce inelastic losses by varying the fre- 
quency of the tight confinement opens new handles on 
tunable BEC in quasi2D gases. These experiments can 
be combined with measurements of non-trivial phase co- 
herence properties of condensates with fluctuating phase. 
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